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Abstract. The paper is devoted to the analysis of the systematic deviations between the classical formula for the 
rate of secular apsidal motion in close binaries and the formula established within the framework of the theory of 
dynamic tides. The basic behaviour is found to be governed by the dependency of the forcing frequencies on the 
orbital period and the rotational angular velocity, irrespective of the stellar model. For binaries with short orbital 
periods a good agreement between the classical apsidal-motion formula and the formula established within the 
framework of the theory of dynamic tides is only possible for certain values of the rotational angular velocity. 
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1. Introduction 

In close binary systems of stars, the perturbation of the 
external gravitational held due to the tidal and rotational 
distortions of the component stars yields a periodic shift of 
the positio n of the periastron. The effect was hrst st udied 
by Russell (1928), and subsequently by Cowling ( 1938 ) 
and ^ter'iiq ( |1939[) . The latter author derived a formula 
for the rate of secular apsidal motion based on the as¬ 
sumption that the orbital period is long in comparison to 
the periods of the free oscillation modes of the component 
stars. The validity of this assumption has been investi- 
gated extensively by, e.g., Papaloizou fc PringL 


Smeyers et al. (1991), iQuataert et al.| (|199( 

3), pmeyers 

et al. (1998|), and Smeyers & Willemg| (^001). 

Smeyers & 


WillemsI (pOOl) found the classical apsidal-motion formula 


One of the most uncertain parameters in the analysis 
of the dynamics of a close binary is the rotational angular 


velocity of the component stars. Claret & Willems (2002) 


to be valid up to high orbital eccentricities as long as the 
orbital and rotational angular velocities remain small. For 
binaries with higher orbital and rotational angular veloc¬ 
ities, the apsidal-motion rates predicted by the classical 
formula deviate from those predicted by the theory of dy¬ 
namic tides due the effects of the compressibility of the 
stellar fluid and the possibility of resonances between dy¬ 
namic tides and free oscillation modes of the component 
stars. 


therefore considered the effects of dynamic tides on the 
apsidal-motion rate as a function of the rotational angu¬ 
lar velocity for a carefully selected sample of eclipsing bi¬ 
naries with accurate determinations of the stellar masses 
and radii and the apsidal-motion rate. The authors found 
that even for binaries with short orbital periods, a good 
agreement between the classical apsidal-motion formula 
and the formula established within the framework of the 
theory of dynamic tides is possible for certain values of 
the rotational angular velocity. 

In this research note, we aim to shed more light on 

(2002) by exam- 


the results found by Claret & Willem: 
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ining more closely the role of the forcing frequencies and 
the effects of the compressibility of the stellar fluid on 
the rate of secular apsidal motion due to the tidal distor¬ 
tions of close binary components. In Sects. 2 and 3, we 
recall the basic assumptions and the equations governing 
the rate of secular change of the longitude of the peri¬ 
astron in the classical framework and in the framework 
of the theory of dynamic tides. In Sect. 4, we illustrate 
the behaviour of the relative deviations between the clas¬ 
sical apsidal-motion formula and the formula established 
within the framework of the theory of dynamic tides as a 
function of the rotational angular velocity in the case of 
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an evolved SMq main-sequence star. In Sects. 5 and 6, 
we unravel this behaviour by looking at the effects of the 
forcing frequencies and the stellar model separately. In the 
final section, we briefly summarise our conclusions. 


2. Basic assumptions 

Consider a close binary system of stars with masses Mi 
and M 2 that are orbiting around each other under the 
influence of their mutual gravitational force. Let Porb be 
the orbital period, a the semi-major axis, and e the orbital 
eccentricity. We assume the first star to rotate uniformly 
around an axis perpendicular to the orbital plane with an 
angular velocity H whose magnitude is small in compari¬ 
son to the stellar break-up speed He- The companion star 
is treated as a point mass. 

Furthermore, we denote by r = (r, d, (f) a system of 
spherical coordinates with respect to an orthogonal frame 
of reference that is corotating with the star. The poten¬ 
tial governing the tidal force exerted by the companion 
then be expanded in terms of 1,1 nnormalised sr 


where w is the longitude of the periastron, ^2 is the clas¬ 
sical apsidal-motion constant, and 
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7 : e 


( 3 ) 


A formula for the rate of secular change of the lon¬ 
gitude of the periastron taking in to account the effects 
of dynamic tides was derived by Smeyers et al.| ( |1998 ) 
by adding the contributions to the apsidal-motion rate 
stemming from the various terms in Expansion (|^ of the 
tide-generating potential. The resulting formula takes the 
form 
dxu\ 

dt )dyn~ 


a 


M 2 2 TT 


Ml Porb 


2 k2 G2 ,o,o + 4 ^ 

fe=l 


(P2,0,fe G2,0,k + P2,2,fe G2,2,k + F2.-2,k G2,-2,k) 


( 4 ) 


spheri- 


can 

cal hailmonics and in Fourier series in terms of 

multiples of the companion’s mean motion n as 

erW{r,t) = -sr^ti: E 


where the constants P 2 ,m,fc render the response of the star 
to the various forcing angular frequencies ur, and the co¬ 
efficients G 2 ,m,k are functions of the orbital eccentricity. 
For a definition of the constants P 2 ,m,fe and the coefficients 
G 2 ,m.k, we refe r to Bmeyers et al. ( 1998| ) and Smeyers & 
Willems! (^00 ij). The coefficients G 2 ,m,k decrease rapidly 


Pi 


£=2 m——i k— — c 


ci,m,k 1 ^ ) Yp{0,(j)) exp[i{aTt - kur)] 


( 1 ) 


(e.g. Polfliet fc Smeyers 199C , ^meyers et alj 199^ ). In this 

expa.ns inn, H is the Newtonian constant of gravitation , /? 


with increasing values of k. 

For binaries with short orbital periods or rapidly rotat¬ 
ing component stars, the rates of secular apsidal motion 
predicted by the classical formula deviate from the corre¬ 
sponding rates predicted by the formula established in the 
framework of the theory of dynamic tides due to the in¬ 
creased role of the stellar compressibility at higher forcing 
frequencies and due to resonances of dynamic tides with 
free oscillation modes of the component stars (Smeyers & 


is the iqnilibriiim ra.diiis of the uniformly rotating star, Willems 2001 ). The extent of the deviations depends on 


Et = {Ri/ay M 2 /M 1 is a small dimensionless parameter, 
r is a time of periastron passage, and ar = kn + mfl is 
a forcing angular frequency with respect to the corotating 
frame of reference. For a definition of the Fourier coeffi- 


the evolutionary stage of the star due to the increase of the 
stellar radius and the redistribution of the mass as the star 


cients ce^m,k, we refer to jSmeyers et al.| ( 11998]) . They are 
equal to zero for odd values of i |m| and, for a given 
eccentricity, decrease rapidly with increasing values of k. 

For the remainder of the paper, we restrict ourselves to 
the contributions of the dominant terms associated with 
^ = 2 in the expansion of the tide-generating potential. 


3. The apsidal-motion rate due to the tidal 
deformations of binary components 


evolves on the main sequence (Willems & Claret 2002). In 
the present investigation, we focus on the role of the forc¬ 
ing frequencies and the compressibility of the stellar fluid. 


In analogy to previous investigations (Smeyers & Willemi 
^001 , Willems & Claret 2002, Claret & Willems 2002)^ 


we study the deviations between the apsidal-motion rates 
predicted by the classical formula and the rates predicted 
by the formula taking into account the effects of dynamic 
tides by means of the relative difference 

classical “ 


A = 


[dw/dt)^^^ 


( 5 ) 


The most commonly used formula for the rate of secu¬ 
lar change of the longitude of the periastron due to the 
tidal deformati ons of the com pone nts of a close binary 
was derived by [Cowling ( 1938 ) and ^terne ( 1939| ) on the 
assumption that the orbital period is long in comparison 
with the free harmonic periods of the component stars. 
The rate of secular apsidal motion is then given by 


dw 

dt 


classical 


l Ri\ ^2 2 7r , 2\ (r)\ 

= (t) mTTuAC. (2) 


4. Relative deviations for a SMq main-sequence 
star 

We illustrate the behaviour of the relative differences A 
as a function of the rotational angular velocity 12 in the 
case of a 5 Mq main-sequence star with a central hydrogen 
abundance Xc = 0.36 and a radius Ri = 3.8 Rq. Since the 
rotational angular velocity 12 is assumed to be small, we 
neglect the effects of the Coriolis force and the centrifugal 
force for the determination of the constants F 2 ^m,k- The 
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Fig. 1. The relative differences A for the evolved 5 Mq 
main-sequence model, the orbital eccentricity e = 0.4, and 
the orbital periods Porb = 5-0, 7.5,10.0,12.5 days. 


tides then correspond to those of a non-rotating spher¬ 
ically symmetric equilibrium star (for more details see, 
e.g., [Willems fc Claret 2002| ). 

The resulting relative deviations A are presented in 
Fig. 0 for the orbital eccentricity e = 0.4 and for orbital 
periods ranging from 5.0 to 12.5 days. For convenience, 
the rotational angular velocity is expressed in units of 
the orbital angular velocity flp at the periastron of the 
relative orbit. The peaks in the curves correspond to reso¬ 
nances between dynamic tides and free oscillation modes. 
Since we are here interested in the systematic trend caused 
by the effects of the forcing frequencies and the compress¬ 
ibility of the stellar fluid rather than in a detailed treat¬ 
ment of resonant dynamic tides, we did not attempt to 
find all possible resonances for each orbital configuration. 
For a more detailed treatment of the resonances, we refer 


the reader to Smeyers & Willems (2001) and Willems & 
Claretf] ([20021) . 

In the limiting case where 11 —> 0, the forcing frequen¬ 
cies ctt are determined solely by the mean motion n. As 
shown by Smeyers & Willems (2001), the relative differ¬ 
ences A are then generally negative and larger in absolute 
value for shorter orbital periods. For larger values of H, the 
relative differences A increase, reach a positive maximum, 
and subsequently decrease to become negative. Hence, we 
find two values of the rotational angular velocity for which 
the relative differences A become zero. These values de¬ 
pend on both the period and the eccentricity of the orbit, 
and do not necessarily correspond to synchronised rota¬ 
tion with the orbital angular velocity at the periastron of 
the relative orbit. We will further scrutinise this behaviour 
in the following section. 


5. The role of the forcing frequencies 


Due to the expansion of the tide-generating potential in 
Fourier series in terms of multiples of the companion’s 
mean motion, the tidally distorted star is subjected to 
an infinite number of forcing angular frequencies ut [see 
Eq. (^]. For a given orbital eccentricity, the largest dif¬ 
ferences between the classical apsidal-motion formula and 
the formula established within the framework of the the¬ 
ory of dynamic tides can be expected to occur when high 
forcing frequencies are associated with large Fourier coef¬ 
ficients C 2 ,m,fe- Vice versa, the smallest differences can be 
expected to occur when the dominant terms in the expan¬ 
sion of the tide-generating potential all have small forcing 
frequencies. 


A similar reasoning applies to Expansion ([j) which 
renders the apsidal-motion rate within the framework of 
the theory of dynamic tides. Since the system of differ¬ 
ential equations governing non-resonant dynamic tides in 
components of close binaries depends on the square of the 
forcing frequency, the deviations of the constants F 2 ^m,k 
from the classical apsidal-mo tion constant fc 2 are of the 
ord er of cr^. Because of thi s, Smeyers fc Willems ( 2001 ) 
and Claret fc Willems ( 2002 ) were able to approximate the 
systematic deviations between the classical apsidal-motion 
formula and the formula established within the framework 
of the theory of dynamic tides by second-degree polyno¬ 
mials in the companion’s mean motion n or the star’s ro¬ 
tational angular velocity 12. In the present context, this 
property implies that the differences between the two for¬ 
mulae can be expected to be large when high forcing fre¬ 
quencies are associated with large coefficients G 2 ,m,k- 


In order to shed some more light on this, we examine 
the behaviour of the function 


2 oo 

5l(Porb,12)= ^ (6) 

m—— 2 

The function mimics the behaviour described above in the 
sense that it will be large when high forcing frequencies 
are associated with large coefficients G 2 ,m,k- 

The variations of Si as a function of the rotational an¬ 
gular velocity 12 are displayed in Fig. in the case of the 
orbital eccentricity e = 0.4 and orbital periods ranging 
from 5.0 to 12.5 days. For given values of e and Porb, the 
function 5'i initially decreases with increasing values of 12. 
The decrease is caused by the decrease of the forcing fre¬ 
quencies associated with m = — 2, which generally provide 
the dominant contributions to 5'i. As 12 becomes larger. 
Si reaches a minimum when the forcing frequencies associ¬ 
ated with the largest coefficients G 2 - 2 ,k are close to zero. 
The frequencies associated with m = —2 subsequently be¬ 
come more and more negative, resulting in an indefinite 
increase of the function 5'i. 

As a preliminary conclusion we can thus say that the 
opposite of the function shows the same qualitative 
behaviour as the relative differences A shown in Fig. |]. In 
the next section, we will make this even clearer by adding 
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n/r}p n/Qp 


Fig. 2. The variation of Si (Porb, as a function of the 
rotational angular velocity for the same orbital eccen¬ 
tricity and orbital periods as considered in Fig. 


Fig. 3. The variation of S 2 (Porb, for the evolved 5 Mq 
main-sequence model and the same orbital eccentricity 
and orbital periods as considered in Fig. |^. 


some information on the stellar model to the definition of 
Pi. 


6. The role of the stellar model 


The response of a particular stellar model to the tidal forc¬ 
ing of a close companion is determined by the constants 
P 2 ,m,fc- For low-frequency, non-resonant dynamic tides in 
stars with a radiative envelope, these constants can be 
approximated by their asymptotic representation 


pi®-) _ i, 

^2,m,k - «2 - 


Pf CTy 
GMi T2 


m 


+ 2 


Ri 


6 (Pi) 
Pi 


(7) 


where ^2 {r) is a solution of the homogeneous second-order 
differential equation 


d'^^2 

dr^ 


1 dm{r) 
m(r) dr 


^-4^2 = 0 (8) 

( ^meyers et al. |l998 ). In the latter equation, m(r) is the 
mass contained inside the sphere with radius r. The solu¬ 
tion of Eq. (|^) must furthermore remain finite at r = 0 
and satisfy the boundary condition 


dr 


+ 


Ri 


6 (Pi) 
Pi 


= 5. 


(9) 


The asymptotic approximation for the constants 
P 2 ,m,fc is valid for forcing frequencies ar ^ \/GMi/Rf. 
Higher forcing frequencies, such as those associated with 
large values of k, require additional terms in the asymp¬ 
totic approximation. However, the second-order approx¬ 
imation given by Eq. (0) is sufficient to further clarify 
the qualitative behaviour of the relative differences A ob¬ 
served in Fig. ^ 


We proceed in a similar way as for the introduction of 
the function Pi in the previous section, but now take into 
account the structure of the stellar model by defining a 
function S 2 as 

2 00 

P2(Porb,f2)= ^ (10) 

m= — 2 k—0 

The function can be rewritten in terms of Pi as 

2 c» 

P2(Porb,F!)=fc2 ^ 

m——2 k—0 


1 pf 

12 GMi 


I g (Pi) 

dr ) „ Pi 


Pi (Porb,f2). 


( 11 ) 


For a given stellar model and a given orbital eccentricity, 
the first term is then seen to be constant, while the sec¬ 
ond term depends on the orbital period and the rotational 
angular velocity through the function Pi. 

The variations of P 2 as a function of the rotational an¬ 
gular velocity Vt are displayed in Fig. ^ for the same 5 Mq 
main-sequence model used in Sect. 4. The qualitative be¬ 
haviour of P 2 is essentially the same as that of —Pi. In 
addition, the inclusion of the model dependent factors in¬ 
troduces two zeros in the curves, similar to the behaviour 
of the relative differences A displayed in Fig. |^. 


7. Summary 

In this research note, we scrutinised the behaviour of the 
systematic relative deviations between the classical for¬ 
mula for the rate of secular apsidal motion and the for¬ 
mula established within the framework of the theory of 
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dynamic tides as a function of the often uncertain rota¬ 
tional angular velocity of the component stars. We found 
that when the effects of resonant dynamic tides are ne¬ 
glected, the basic behaviour of the relative differences is 
governed by the dependency of the forcing frequencies on 
the orbital period and the rotational angular velocity, ir¬ 
respective of the stellar model. For a given value of the 
orbital eccentricity, a good agreement between the clas¬ 
sical apsidal-motion formula and the formula established 
within the framework of the theory of dynamic tides is 
possible even for binaries with short orbital periods, but 
only for certain values of the rotational angular velocity. 
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